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ABSTRACT 

Using the results for the NNLO QCD corrections to anomalous dimensions of odd xF^ Mellin moments and N^LO 
corrections to their coefficient functions we improve our previous analysis of the CCFR'97 data for xF^. The possi- 
bility of extracting from the fits of 1/Q^-corrections is analysed using three independent models, including infrared 
renormalon one. Theoretical question of applicability of the renormalon-type inspired large-/3o approximation for 
estimating corrections to the coefficient functions of odd xF^ and even non-singlet F2 moments are considered. The 
comparison with [1/1] Pade estimates is given. The obtained NLO and NNLO values of as{Mz) are supporting the 
results of our less definite previous analysis and are in agreement with the world average value as{Mz) ~ 0.118. We 
also present firts N^LO extraction of as{Mz). The interplay between higher-order perturbative QCD corrections and 
l/Q^-terms is demonstrated. The results of our studies are compared with those obtained recently using the NNLO 
model of the kernel of the DGLAP equation and with the results of the NNLO fits to CCFR'97 xF^ data, performed 
by the Bernstein polynomial technique. 
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1 Introduction 



The series of our previous works [l]-[5] (for a brief summary see the review of Ref. [6]) was devoted to the QCD 
analysis of the xF^ data obtained at the Fermilab Tevatron. In Refs. [3]- [5] we made the consequent steps towards 
the next-to-next-to-leading order (NNLO) fits to the xF^ CCFR'97 data [7] both without and with twist-4 corrections 
taken into account. In the process of these studies definite theoretical approximations were made. This work is devoted 
to the fixation of the number of theoretical ambiguities involved in our previous fits. Let us recall the basic steps of 
these studies. 

The a;-behaviour of xF^ was reconstructed from the Mellin moments with the number 2 < n < 10 with the help 
of the Jacobi polynomial technique, developed in Refs. [8] -[12]. At the NNLO the perturbative expansion for the 
coefficient functions of these moments is explicitly known from analytical calculations of Refs. [13], which were con- 
firmed recently with the help of another technique [14]. However, the NNLO corrections to the anomalous dimensions 
of the considered non-singlet (NS) moments, closely related to still analytically unknown NNLO contributions to the 
kernel of the DGLAP equation [15]-[18], were modelled in the definite approximations only. These approximations 
were based on the observation that for 2 < n < 10 the next-to-lcading (NLO) corrections to the anomalous dimensions 
f-2 ('^) ^^"^ Mellin moments for F2 do not numerically differ from the NLO corrections to the anomalous 
dimensions 7|r^''(?T-) of the moments taken from the xF^ structure function (SF) [1]. In view of this it was assumed 
that this similarity will be true in the case of higher-order corrections to the anomalous dimensions and 7^^^^ 

as well, provided the typical diagrams with new Casimir structure d"^'^d°'^° (which are starting to contribute to 7^^ 
from the NNLO) are not large. Using this assumption, we applied in our analysis [3]- [5] of the xF^ CCFR'97 data [7] 
the available results for the NNLO corrections to 7]vsf2('^)' which were known in the cases of n = 2,4,6,8, 10 due 
to the distinguished calculations of Refs. [19, 20]. For the odd moments with n = 3,5, 7, 9 the NNLO corrections to 
INS F- (^*) were estimated using the procedure of smooth interpolations, previously proposed in Ref. [21] in the process 
of first NS NNLO fits to the F2 data of the BCDMS collaboration [22]. 

Quite recently the renormalization group quantities for the xF^, Mellin moments were analytically calculated at 
the NNLO level [23] by the methods of Refs. [24]- [26] , used in the case of calculations of the even Mellin moments of 
F2 [19, 20]. In Ref. [23] the following information, quite useful for the fixation of some theoretical ambiguities of our 
previous fits, was obtained: 

in) 

• the NNLO corrections to ^^p' at n = 3, 5, 7, 9, 11,13; 

in) 

• the NNLO corrections to Tns,F2 at n = 12 and n = 14; 

• the N^LO corrections to the coefficient functions of odd Mellin moments of xF^ with n = 3, 5, 7, 9, 11, 13 ^. 

Due to the appearance of this new important information it became possible to recast the analysis of Refs. [3]- [5] on 
the new level of understanding. First, we are now able to improve the precision of the smooth interpolation procedure 
for J^j^g p^; applied in the process of the works of Refs. [l]-[5]. Second, it is now possible to change the previously- used 
approximation for NS anomalous dimensions and take into account the available expressions for the NNLO corrections 
to "fFsin-) including terms typical to this level of perturbation theory, which are proportional to /(c?"'"^(i"^^)/n (in our 
case we will consider / = 4 number of massless active flavours, while for the 5C/c'(3)-group ((i"''^(i"''^)/3 = 40/9). 
Moreover, since the best way of estimating the perturbative part of theoretical uncertainties of the N^^^LO fits is the 
incorporation of the N^'+^^LO-terms, the results of Ref. [23] are giving us the chance to perform the approximate 
N'^LO analysis of xF^ data and compare its outcome with the results obtained within the framework of the N^LO 
Fade motivated fits of Ref. [5] . 

The aims of this work are the following: 

1. To reveal and eliminate definite theoretical uncertainties related to the previously- used NNLO approximation 
7^V2 W ~ ^F? W at 2 < n < 10; 

2. To include in the analysis the NNLO corrections to the coefficient functions Cp^ and to the anomalous dimensions 
7j^^ up to n < 13 (the first results of this program were already presented in Ref. [28]); 

^In the case of n = 1 the contribution to the Gross-Llewellyn-Smith sum rule was already known [27]. 



3. To reveal and fix the uncertainties of our previous NNLO fits of Refs. [3]- [5], which had the aim to determine 
the NNLO value of as{Mz) and the parameters of the xF^ SF model at the initial scale Q^; 

4. To check the feature of the interplay between NNLO perturbative QCD corrections and terms, discovered 
in Refs. [3, 5], which leads to the cflFcctivc decrease of the value of the basic free parameter of the infrared 
renormalon (IRR) model for the twist-4 contribution to xF^ [29] (for the discussions of the IRR contributions 
to characteristics of deep-inelastic scattering see Refs. [30, 31], while for recent reviews of the current general 
status of the IRR approach see Refs. [32, 33]); 

5. To comment upon the predictive abilities and special features of the Fade approximants at the N^LO level; 

6. To study the applicability of the renormalon-inspired large- /3o expansion for estimating higher-order perturbative 
corrections to the coefficient functions of moments of xF^, which will be used by us. The similar consideration 
of NS moments for F2 SF [34] is also updated. The results will be compared with those obtained by means of 
the Fade approximation technique; 

7. To study the scale-dependence of the obtained values for Ug at the NLO and beyond; 

8. To reconsider the problem of the extraction of the a;-shape of the twist-4 contributions at the LO, NLO, NNLO 

(2) 

using new information about 'ypj{n) and taking into account additional terms in the perturbative expansion of 
dpj at the N^LO; 

9. To compare our results with other phenomenological applications of the perturbative QCD calculations of 
Refs. [13, 23], which appeared in the literature recently and were based on the application of DGLAP approach 
(see Ref. [35]) and Bernstein polynomial technique proposed in Ref. [36] and used at the NNLO in Ref. [37]. 



It should be stressed that since in the current work we are interested in the study of theoretical ambiguities of 
the outcome of the fits on a more solid background than was done in Ref. [5], and in revealing the importance of 
the knowledge of exact expressions for still uncalculated NNLO corrections to Jpj{n) for n even, we will neglect 
the systematic experimental uncertainties of the CCFR'97 xF^ data, taking into account statistical ones only. The 
incorporation to the fits of systematic error-bars might shadow the effects of fixing the theoretical uncertainties we 
are struggling for and might make the process of clarification of the necessity of getting explicit numbers for 7^^-' (n) 
for even n more complicated. At the NLO level the combined analysis of statistical and systematic experimental 
uncertainties of the xF^ CCFR'97 data was done in Refs. [38, 39] using the machinery of the method proposed in 
Ref. [40]. 



2 Preliminaries 

For the sake of completeness of the presentation we will repeat some definitions from Ref. [5] . 
We start from the Mellin moments for 2:^3(2;, Q'^): 

M^^{Q^)= [\"-^F3{x,Q^)dx (1) 
Jo 

where n = 2, 3, 4, .... These moments obey the following renormalization group equation 

+/3(^)^ -7l.:^(^))Ae(QVM',^(M')) = (2) 
where Ag = as/(47r). The renormalization group functions are defined as 

/i^=/3(^) = -25^/3,4+2 

i>0 



-M^=7S^(A) = E7g(nK^^ (3) 

" i>0 



The NS anomalous dimensions of F2 will be defined in the analogous way, namely : 



i>0 



where Z^^ and Z^^--^^ are the renormahzation constants of the corresponding NS operators. In the case of xF^ 
moments the solution of the renormahzation group equation is: 



exp 



(5) 



(6) 



This expression is identical to the parametrization used by the CCFR collaboration [7] . In the process of our studies 
we will extract from the fits the parameters A{Qq), 6(Qo)) c((5o) and 7((5o)) together with the parameter A^IL and 
with the information on the twist-4 terms. 

The first moment of xFs coincides with the Gross-Llewellyn Smith sum rule 



where M^^{Qq) is a phenomenological quantity defined at the initial scale Qq as: 

M^'{Ql) = C x''-^A{Ql)x'^^o)^i-:,yiQl)(i + ^{^Ql):,)dx . 



GLS{Q^) = M['{Q^) = t Faix, Q^)da 

Jo 



At the N'^LO the coefficient function Cj^'', which enters Eq. (5), can be defined as 

C^J{As) = 1 + cg\n)As + C^^{n)Al + C^^{n)Al 
while the corresponding expression for the anomalous-dimension term is : 



exp 



p{x) 



-dx 



{A,{Q')y^>^/'^°xAD{n,A,) 



where 



AD{n,A,) = [l+p{n)A,{Q^) + q{n)A,{Q^f+r{n)A,{Q^f] 
and p{n), q{n) and r{n) read: 



P{n) = - 



/3i 



/3o ;/3o 



q{n) 



2p{nf + 



r(n) = -^l - 2p{nf + 6p{n)q{n) + 



7g(n) /3i7g(n) 



f3o 



Pi 



(7) 

(8) 

(9) 
(10) 

(11) 
(12) 

(13) 



P2li\n)^Ph^^{n) Ph^Si^) PilfM 2/31/3271.? (n) 



PI 



P'o 



Pt 



0i 



PI 



The coupling constant ^s(Q^) can be decomposed into the inverse powers of i = ln((5^/A|^) as A^^'~' = A^*^ + 
A^f ^o, A^^^o = j^NLO + AA^^^o and Af^^ = j^nnlo + ^Af^o^ ^i^g^e 

A^o 

1 



1 

PilnjL) 
PlL^ 



NNLO 



PIL^ 



[Piln'iL) - /3tln{L) + P2P0 - Pf] 



aa: 



^^[(3l{-ln\L) + ^In^L) + 2ln{L) - \) 

-Zl3oPip2ln{L) + Pl^] . 



(14) 

(15) 
(16) 



Notice that in our normalization the expressions for /Jq, Pi, P2 and have the following numerical expressions: 

/Jo = 11 - 0.6667/ 

/3i = 102 - 12.6667/ 

(32 = 1428.50- 279.611/ + 6.01852/2 

/?3 = 29243.0 -6946.30/ + 405.089/2 + 1.49931/^ (17) 
where P3 was analytically calculated in Ref.[41]. 



3 Anomalous dimensions and coefficient functions: 
approximate vs exact results 

Let us discuss the numerical approximations of higher-order perturbative corrections to anomalous dimensions and 
coefficient functions, used in our fits. The analytical expression for the one-loop term of NS anomalous dimensions 
"fpf^sp^in) = "/p^{n) = 8/3[4^"^j^(l/j) — 2/n(n+ 1) — 3] is well known. The NLO corrections to the NS anomalous 
dimensions were obtained in Refs. [42, 43] and confirmed by the independent calculation in Ref. [44]. In the cases of 
both F2 and xF^, the numerical expressions for the NLO contributions to the anomalous dimensions are given in Table 

1, where we also present the numerical values of the NNLO coefficients used to the NS anomalous dimension functions 

(2) 

of Eq. (3) and Eq. (4). They are also normalized to the case of / = 4 number of flavours. The coefficients Jp,' {n)\wts 
represent the contribution to 7^' of the terms without typical structure (wts) f {d"'''^ d"''"^) / n = 4*40/9, which is absent 
in the expression for 7]^^ and appears for the first time in the anomalous dimensions of XF3 moments at the NNLO. 

We describe now how the related approximations for 7^^ (n) and 7^'' (n) were obtained and what is the connection 
between the sets of these numbers. The expressions without round brackets are the results of explicit analytical 
calculations of Refs. [19, 20, 23]. The expressions in round brackets arc the approximations obtained with the help of 
the smooth interpolation procedure. To study the possibility of getting stable dependence of the values for from 

the change of the initial scale Qq after incorporation to the ffis of the NNLO corrections to 7^''' with 10 < n < 13, the 

(2) 

interpolation procedure was supplemented with the fine tuning of the NNLO corrections to Jpj {n) with n = 6, 8, 10. 



n 








7f, (")Ut« 




2 


71.374 


71.241 


612.0598 


(585) 


(631) 


3 


100.801 


100.782 


( 839.8534) 


836.3440 


861.6526 


4 


120.145 


120.140 


1005.823 


(1001.418) 


(1015.368) 


5 


134.905 


134.903 


(1134.967) 


1132.727 


1140.900 


6 


147.003 


147.002 


1242.0006 


(1241.21) 


(1246.59) 


7 


157.332 


157.332 


(1334.865) 


1334.316 


1338.272 


8 


166.386 


166.386 


1417.451 


(1416.73) 


(1419.783) 


9 


174.468 


174.468 


(1491.711) 


1491.124 


1493.466 


10 


181.781 


181.781 


1559.005 


(1558.854) 


(1560.675) 


11 


188.466 


188.466 


(1620.755) 


1620.727 


1622.283 


12 


194.629 


194.629 


1678.400 


(1677.696) 


(1679.809) 


13 


200.350 


200.350 


(1732.640) 


1731.696 


1732.809 



Table 1. The numerical expressions for NLO and NNLO coefficients of anomalous dimensions of NS moments of F2 
and xFs at / = 4 number of flavours. 

The numerical values of the obtained fine-tuned numbers will be given in the next Section. This procedure results 
in better stability of the fitted values of A^-}— with respect to changes of the initial scale at Qq > 5 GeV^ (see Table 4 

(2) 

below). The approximation for 7)^^ (2) contains more uncertainties. At the first stage it was obtained by extrapolation 

of the NNLO coefiicients to anomalous dimensions with n> 2 without using the explicit number for 7^^ (n = 1) (which 

is zero). We have checked the reliability of this procedure, considering the set of NLO anomalous dimensions jp^{n) 
which are explicitly known at any values of n. We found that the interpolated values at even n = 4,6, ... are much 
closer to the real numbers than those obtained with incorporation of the the zero expression for the n = 1 anomalous 



dimension. In this case for n = 2 we obtained the extrapolated value (75.41), which is over 6% higher than the real 
value 71.24. In the case of application of the interpolation procedure, which used the n = 1 result (zero), we get the 
value (66.11), which is 6% smaller than the real value. To estimate the NNLO expression for ^p^{2) we imposed the 
conditions of the reduction found at the NLO and thus fixed the value of (2) 6% below the extrapolated number. 
Its numerical expression, which will be used throughout this work, is quoted in Table 1. The 6% uncertainty of jpj (2) 

translates to a 1 — 2 MeV variation of A^IL, which is below the precision of the NNLO extraction of the value of this 
parameter. 

In view of the doubts in the validity of the interpolation procedure, expressed in Ref. [45], it is of definite interest 
to desribe the results of its applications in more detail. At the first stage one can compare the interpolated numerical 

expressions for jj^g p^{n) at n = 3, 5, 7, 9, 11, 13, obtained with application of values recently calculated in Ref. [23] 

for 7/^5^ir2(12) and Jn\^F2(^'^) with the explicit identical results for 7^^(n)|u)ts) which do not contain the cl°'''c^abc_ 
structure. The comparison is presented in Table 1. The estimates from the third column (839.8534), (1134.967), 
(1334.865), (1491.711), (1620.755), (1732.640) should be compared with the explicit numbers 836.3440; 1132.727; 
1334.316; 1491.124; 1620.727 and 1731.696. One can see that the qualitative agreement between these sets of numbers 
is rather good. Moreover, we can study the applicability of the interpolation procedure for simulating explicitly 
unknown coefficients of jj^fg FaC"-) obtained with and without application of j}^g p2{^2) and 7]^^ p2{^^) terms. In fact 
we found that for ri=3,5,7,9 the difference between the "new" and "old" interpolated numbers is not large, namely 
2.45; -0.8; 0.9; -1.8. It affects the third significant digit of the interpolated estimates and improves their qualitative 
agreement with the results of the calculations of ^pj{n) at odd values of n [23]. However, it is known that even a 
knowledge of the 3rd significant digit in anomalous dimension terms is not enough for the precise reconstruction of 
SF from the the NLO results for the moments at large n > 6 [11, 12]. Therefore, to determine the unknown even 

f2) (2.) 

values of jj^g p^{n) from the known NNLO corrections to 7)1-3 (n) at n = 3, 5, 7, 9, 11, 13 [23] with more precision we 

supplement the interpolation procedure for ^p^{n) by the fine-tuning of its terms for n = 6,8,10. The results of 

the application of the first approximation are presented in the last column of Table 1 . Notice that they contain the 
scheme-independent new contributions, labelled by the d'^^'^d"'^'^ gauge group structure. The fine-tuned expressions for 
"fpjin) at n = 6, 8, 10 will be presented below. 

One more verification of the idea of smooth interpolation comes from the consideration of its application for 
estimating NNLO corrections to the coefficient functions of even Mellin moments of xF^. The estimates obtained are 
presented in column 4 of Table 2. The comparison with the results, given in column 3, which were calculated from 
the expression of Ref. [13], demonstrate perfect agreement of these estimates with the explicit numbers. We think 
that in view of this, one can safely apply the idea of smooth interpolation in order to estimate the N^LO terms to the 
coefficient functions of even Mellin moments of xF^ from those calculated in Ref. [23] with n = 3, 5, 7, 9, 11, 13, and 
also taking into account the order 0{Al) correction to the Gross-Llewellyn Smith sum rule (n = 1 moment), obtained 
in Ref. [27]. 

The information about the N^LO corrections to the considered coefficient functions is the important and dominating 

ingredient of the N'^LO fits to xF^ data we are going to perform. It should be supplemented with the model for the 
N^LO corrections r(n) to the anomalous dimension (AD) function of Eq. (10). 



n 






Cl:'i.)U. 




^1'' ('0,11/1, 




1 


-4 


-52 


-52 


-644.3464 


-676 


480 


2 


-1.778 


-47.472 


(-46.4295) 


(-1127.454) 


-1267.643 


174.4079 


3 


1.667 


-12.715 


-12.715 


-1013.171 


97.00418 


-47.01328 


4 


4.867 


37.117 


(37.0076) 


(-410.6652) 


283.0851 


246.0090 


5 


7.748 


95.4086 


95.4086 


584.9453 


1174.835 


1013.328 


6 


10.351 


158.2912 


(158.4032) 


(1893.575) 


2420.569 


2167.903 


7 


12.722 


223.8978 


223.8978 


3450.468 


3940.284 


3637.790 


8 


14.900 


290.8840 


(290.8421) 


(5205.389) 


5678.657 


5360.371 


9 


16.915 


358.5874 


358.5874 


7120.985 


7601.721 


7291.305 


10 


18.791 


426.4422 


(426.5512) 


(9170.207) 


9677.391 


9391.308 


11 


20.544 


494.1881 


494.1881 


11332.82 


11885.25 


11633.28 


12 


22.201 


561.5591 


(561.2668) 


(13590.97) 


14204.22 


13991.80 


13 


23.762 


628.4539 


628.4539 


15923.91 


16620.99 


16449.68 



Table 2. The values for NLO, NNLO, N'^LO QCD contributions to the coefficient functions, used in our fits, and 



the results of N^LO Pade estimates. We fix them using the [1/1] Pade resummation procedure of the coefficients of 

the AD-function (for the results see Table 3, whore the expressions which come from [0/2] Pade estimates arc also 
presented). It should be mentioned that the values for p{n) and q{n) are calculated from the numbers given in Table 
1. 



n 


p{n) 


q{n} 




^Hl[0/2l 


1 














2 


1.6462 


4.8121 


14.0666 


11.3822 


3 


1.9402 


5.5018 


15.6011 


14.0456 


4 


2.0504 


5.8327 


16.5919 


15.2986 


5 


2.1149 


6.2836 


18.6691 


17.1187 


6 


2.1650 


6.7445 


21.0110 


19.0560 


7 


2.2098 


7.1671 


23.2447 


20.8847 


8 


2.2525 


7.6013 


25.6518 


22.8152 


9 


2.2939 


8.0164 


28.0151 


24.7073 


10 


2.3344 


8.4353 


30.4804 


26.6614 


11 


2.3741 


8.8146 


32.7261 


28.4720 


12 


2.4131 


9.1855 


34.9647 


30.2794 


13 


2.4512 


9.5620 


37.3002 


32.1491 



Table 3. The values for NLO and NNLO QCD contributions to the expanded anomalous dimension terms used in 
our fits and the N'^LO Pade estimates. Several comments should be made concerning the comparison of the Pade 

estimates technique of the N^^LO corrections Cp^{n) with more definite, to our point of view, results of application 
of the interpolation procedure (see Table 2). One can see that the agreement of [1/1] Pade estimates with the N^LO 
coefficients is good in the case of the Gross-Llewellyn Smith sum rule (this fact was already known from the estimates 
of Ref. [46], which are close to the results of the scheme- invariant approach of Ref. [47]). In the case of n = 2 and 
n > 6 moments the numbers of columns 5 and 6 of Table 2 are also in satisfactory agreement. Indeed, one should keep 
in mind that the difference between the numbers presented in columns 5 and 6 of Table 2 should be divided by the 
factor (4)'^, which comes from our definition of the expansion parameter Ag = as/(47r). Note that starting from n > 6 
the results of the application of [0/2] Pade approximants, which in accordance with the analysis of Ref. [48] reduce 
scale-dependence uncertainties, are even closer to the the estimates which are given by the interpolation procedure 
(for the comparison of the outcome of approximate N'^LO fits to xF^ data, which are based on [1/1] and [0/2] Pade 
approximants, see Ref. [5], while the comparison of the applications of [1/1] and [0/2] Pade approximants within the 
quantum mechanic model was analysed in Ref. [49]). For n = 3, 4 the interpolation method gives completely different 
results. The failure of the application of the Pade estimates approach in these cases might be related to the irregular 
sign structure of the perturbative series under consideration. A similar problem arises in the case of the analysis of 
the perturbative series for QED renormalization group functions (for discussions see Ref. [50]). 

However, in the case of the perturbative series AD for the expanded anomalous dimension term we do not face 
this problem (see Table 3). In view of the absence of other ways of fixation of N^LO coefficients r(n) (the renormalon- 
inspired large-/3o expansion is definitely not working for the anomalous dimensions functions [51, 52]) we will use in 
our fits the [1/1] Pade estimates for r{n). Note in advance that the application of [0/2] Pade resummation to r(n) 
does not influence the outcome of our approximate N^LO fits with defined by its explicit N^LO expression (see 
Eqs.(14)-(16)). It should be stressed that since in the process of these fits wc will use the explicitly calculated coefficient 
functions of the xi^a Mellin moments [23] , the obtained uncertainties will be more definite than those estimated in our 
previous work of Ref. [5] . 



4 Results of the fits without twist-4 terms 

In order to perform the concrete fits to xF^ CCFR'97 data and thus analyze how new theoretical input, described in 
Sec. 3, affects the results previously obtained in Ref. [5], we apply the same theoretical method, based on reconstruction 
of xF^ from its Mellin moments using the Jacobi polynomial expansion [8]- [12]: 

xFi'--{x,Q')=w{a,P){x) Y: e:?.''(x)^c5")(a,/J)M™^3(Q^)+^ (18) 

n=0 j=0 ^ 



where Q"^'^ are the Jacobi polynomials, Cj {a,(i) contain a- and /3-dependent Euler F-functions where a,P are the 
Jacobi polynomial parameters, fixed by the minimization of the error in the reconstruction of the SF, and 'w{a, (3) = 
— xY is the corresponding weight function with a = 0.7 and /3 = 3 chosen following the detailed analysis of 
Ref. [5]. The contributions of the dynamical twist-4 terms are modelled by the Q^-independent function h{x). The 
kinematical power corrections, namely the target mass contributions, are included in the reconstruction formula of 
Eq. (??) up to order 0{Ml^^i/Q'^)-terms: 

Ml^Q') = M„^3(Q2) ^ n{n + l) ^1^,^ . (19) 

Using Eqs. (??),(19) one can conclude that choosing Nmax = 6, as was done in the case of NNLO fits of Rcfs. [3]- [5], 
we are taking into account 2 < n < 10 Mellin moments in Eq. (??). As was emphasized above, definite information on 
the NNLO QCD corrections to the renormalization-group functions of 2 < n < 13 moments of xF^ is now available. 
In view of this we can now increase the number of Nmax from Nmax = 6 to Nmax = 9 and analyze the changes in the 
results of the NNLO fits of Ref. [5] due to application of additional theoretical information. 

We are starting our studies from the case when twist-4 contributions are switched off (namely h{x) =0). In Table 
4 we present the dependence of the extracted values of A^l^ from the variations of Nmax and of the initial scale Qq. 
Looking carefully at Table 4, one can clearly see that the results of the NLO fits are rather stable under the changes 
of Qq. Moreover, they are in agreement with the ones presented in Table 3 of Ref. [5]. However, for Nmax = 10 the 
values of are larger than for the case of Nmax = 9. Moreover, we checked that for Qq=20 GeV^ and Nmax=^^, 
despite the fact that the corresponding value — 334 ± 37 MeV is comparable with the one obtained in the case 
of choosing Nm,ax=^Q (see Table 4), continues to increase (we got x^=88.9/86). Therefore, it might be reasonable 
to stop at Nmax=^ and thus take into account in Eq. (18) 13 moments only. 

At Nmax = 6 new NNLO results agree with the findings of Ref. [5] . They demonstrate the same dependence of A^l^ 

from Qq. Notice that it has the stability plateau starting from Qq = 20 GcV^ only. However, there is the important 
difference between the results of the NNLO fits of the current work and the ones of Ref. [5]. Indeed, at the NNLO 
Table 4 demonstrates the widening of the stability plateau for A^i^ to lower Qq values for 7 < Nmax < 9 and the 
minimization of x^-value at Nmax = 9. This welcome feature is showing us the importance of changing the approximate 

expressions for 7]vs(^) from the model ^p.'{n) « y^^'g jr.^{n) used in Ref. [5] to the new one, which is based on the exact 

(2) 

numbers for •-fpj (n), calculated in Ref. [23] plus application of the fine-tuning procedure, which in the case of Nmax ~ 6 
gives "f''p,j{6) = 1247.4222 ±2.1357; for Nmax = 7 gives 7 (6) = 1248.1219 ± 1.0359, 7£'(8) = 1420.1729 ± 4.0854; for 
Nmax = 8 gives j'p^{6) = 1248.5610 ± 1.2951, 7^,3^(8) = 1419.3301 ± 1.5112; j'p^iW) = 1561.4299 ± 1.4074; and for 

Nmax = 9 gives ^'^p^ie) = 1247.7852 ± 0.5091, -f'^p^iS) = 1420.2215 ± 0.3337; 7^3^X0) = 1560.8461 ± 0.2292. Within 

the quoted error bars these fine-tuned numbers agree with the estimates obtained by smooth polynomial interpolation 
and presented in the last column of Table 1. The agreement improves with the increasing of Nmax- In the case of 
Nmax = 9 the difference is in the 4th significant digit. 
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Table 4. The Qq and A^max dependence of aI^L. [MeV]. The values of are presented in parenthesis. 

However, the improved Qo-independent NNLO values of A^^L do not differ significantly from the results of Ref. [5] 
(the difference of over 4-5 MeV is about 7 times smaller than the existing statistical error). A similar feature 
reveals itself in the process of the N^LO fits, which are based in part on application of the exact numbers for the 
cig-corrections to the coefficient functions of odd Mellin moments, calculated in Ref. [23], and were performed using 
the N^LO approximation of ag. It should be stressed that the inclusion in the fits of these numbers at Nmax = 6 does 
not lead to the detectable difference of the new results from the ones obtained in Ref. [5] with the help of the Fade 
approximation method. The essential advantage of the new considerations is that we are able to reach Nmax = 9 and 
observe perfect stability of both NNLO and new N^LO values for A^ to the variation of Qq after slight modification 

of (6), (8), and 7^^(10) (which arise from the application of fine-tuning procedure at the N^LO), and can 
determine from these numbers the Fade approximations for r(n)|[i/i] at n = 6,8,10. It should be mentioned that 

the obtained estimates for the fitted three terms of 7];.^ (n) are in agreement with the estimates presented above at 
the level of 4 significant digits. Note also, that if we use in the fits [0/2] Fade approximants for modelling the N-^LO 
correction r(n), we get the values of A^l^, which are rather close to those obtained in the case of application of [1/1] 
Fade approximants. 

Comparing now the results of the NNLO and approximate N^LO fits we conclude that for 7 < N^ax £ 9 the 
difference between the obatined values of A^IL is rather small and almost disappears for = ^ GeV^. Thus we 
can make the conclusion that we observe the minimization of theoretical uncertainties and, probably, the saturation 
of the predictive power of the corresponding perturbative series at the 4-loop level. A similar feature was discovered 
in th(^ proc;ess of calciilations of the perturbative corrections to the correlator of scalar quark currents in the large-/ 
approximation [53]. Therefore, we think that the perturbation theory approximants for xF;i, moments can be safely 
truncated at one more step beyond the NNLO. Higher-order calculations might manifest the signal for asymptotic 
divergence of the related perturbative QCD predictions and as a result, the increase of the value of x^- 
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Table 5. The determined values of the parameters A, b, c, 7 of the model for xFs and their comparison with the values 
obtained in Ref. [4]. The new ones, related to the NNLO, are marked by bold type. 

Another consequence of our new improved analysis corresponds to the determination of QQ-dependence of the 
parameters A, b, c, and 7. It is presented in Table 5 and is compared with the previous extraction of their Qq- 
depcndence given in Ref. [4]. The LO and NLO numbers are the same, while at the NNLO the new results marked out 
by bold type, which correspond to Nmax = 9, are compared with the previous NNLO ones [4], obtained at Nmax = 6 
using the calculations of 7^^ j;.^ (n)-terms [19, 20] available at this time. One can see the noticeable decrease in at 
the NNLO. In the LO the obtained QQ-dcpcndence of c is in agreement with the shape of variation of this parameter, 
predicted in Ref. [54] and confirmed recently in Ref. [55] . 

Even more interesting to study our NNLO results for the parameter &(Qo)- First, it is almost Qg'iiidependent 
within statistical errors. This fact is in agreement with theoretical demonstration of its Qg'iiidependence, presented 
in Ref. [56] using DGLAP equation. Another result was obtained recently in Ref. [57] in all-orders of 1-loop expression 
for as, using in part the approach, developed in Ref. [58]. It should be stressed that the estimate uj~ = 0.4 in the 

expression Fj^g = (l/x)^ (^)'^ obtained in Ref. [57] for ji « 5.5 GeV, Aqcd = 0.1 and / = 3 is in good numbers 
for 1 — 5 especially in the region > jj,^ = 30 GeV^, which is not considerebly affected by the transformation from 
A(3) » 0.1 GeV, used in Ref. [57], to A^ w 0.4 GeV, as advocated by us. 



5 Incorporation of the twist-4 terms 
5.1 Infared renormalon model pcirameterization 

The next stage in modification of QCD thcorcrctical approximations is the inclusion of the higher-twist terms in the 
expression for the structure functions. At the first stage we will rely on the prediction of the IRR approach [29] and 
model the twist-4 contribution to h{x) in Eq.(18) as 

^ n=0 j=0 

where 

W") = C{n)M^'{Q^)^ with C{n) = -n-4 + 2/(n + 1) + 4/(rz + 2) + 4S,{n). (21) 



The results of the new improved fits to CCFR'97 data for xF^ with the twist-4 term taken into account through 
Eq. (21) are presented in Table 6. 



ordcr/7V„ 



MS 



Gc 



20 GcV^ 



100 GcV^ 



LO/6 
LO/9 



A. 



X /nep 
A' 



A. 



(4^ 
'MS 



X /nep 
A' 



433±54 

81.2/86 
-0.331±0.057 

447±54 

79.8/86 
-0.340±0.059 



431±36 

81.2/86 
-0.330±0.059 

443±54 

80.1/86 
-0.337±0.059 



429±35 
80.6 
-0.328±0.058 

439±56 

79.6/86 
-0.335±0.059 



FT 

MS 



NLO/6 

NLO/9 



A 



X /nep 

A'2 
AllL 

MS 

X^/nep 
A' 



370±38 

80.2/86 
-0.121±0.052 
379±41 
78.6/86 
-0.125±0.053 



369±41 

80.4/86 
-0.121±0.053 

376±39 

79.5/86 
-0.125±0.053 



367±38 

79.9/86 
-0.120±0.052 

374±42 
79.0 
-0.124±0.053 



FT 

'MS 



NNLO/6 
NNLO/7 
NNLO/8 
NNLO/9 



A 



X /nep 



A 



A2 

(4) 



MS 



X /nep 
A'2 

(4) 
MS 



A 



X /nep 

A'2 
AllL 

MS 

X^/ncp 
A' 



297±30 
77.9/86 

-0.007±0.051 
327±34 
75.8/86 

-0.011±0.051 
335±37 
73.8/86 

-0.012±0.051 
331±33 
73.1/86 

-0.013±0.051 



328±36 
76.8/86 

-0.017±0.051 
327±35 
78.1/86 

-0.013±0.051 
330±36 
77.0/86 

-0.013±0.051 
332±35 
75.7/86 

-0.015±0.051 



328±35 
79.5/86 

-0.015±0.051 
328±36 
78.6/86 

-0.015±0.051 
329±36 
77.5/86 

-0.015±0.051 
331±35 
76.9/86 

-0.016±0.051 



'MS 



N3LO/6 
N3LO/7 
N^LO/S 
N3LO/9 



A 



X /nep 

A'2 

(4) 
'MS 



A 



X /nep 
A'2 

MS 

X^/ncp 

a; 



A 



(4^ 
'MS 



X /nep 



305±29 
76.0/86 

0.036±0.051 
331±33 
75.6/86 

0.040±0.052 
337±34 
74.1/86 

0.040±0.052 
333±34 
73.8/86 

0.038±0.052 



327±34 
76.2/86 

0.033±0.052 
325±34 
77.7/86 

0.036±0.052 
326±34 
76.9/86 

0.036±0.052 
328±33 
75.9/86 

0.035±0.052 



326±34 
78.5/86 

0.029±0.052 
326±34 
77.8/86 

0.035±0.052 
326±34 
76.7/86 

0.035±0.052 
328±38 
76.4/86 

0.034±0.052 



Table 6. The results of the fits to the CCFR'97 xF^ data with HT terms modelled through the IRR model. A2' is 
the additional parameter of the fit. The cases of different Qq and N^ax are considered. 



Looking carefully at Table 4 we arrive at the following new conclusions: 



• The x^-value decreases from LO up to NNLO and at N^LO level it almost coincides with the one obtained 

at the NNLO. Moreover, decreases with the increasing of N^ax and distinguishes the fits with Nmax = 9. 
This is the welcome feature of including in the fitting procedure more detailed information on the perturbative 
theory contributions both to coefficient functions and anomalous dimensions of 0:^3 moments, and in partic- 
ular explicitly-calculated three-loop coefficients 7^^(n) at n = 3,5,7,9,11,13 [23], supplemented by us with 
application of the interpolation procedure for even n plus fine-tuning of the terms 7^^(6), (8) and 7^^(10). 

• At the N'^LO, smaller than the one obtained in the process of pure Padc-motivatcd fits of Ref. [5]. This 
feature is related to the fact that the explicit expressions for N^LO corrections for odd XF3 moments [23] are 
now taken into account. 

• For Nmax = 9 the values of and the IRR model free parameter A2 are rather stable to variation of Qq not 
only at the LO, NLO but at the NNLO and N^LO as well. The last property gives favour to our new results in 
comparison with the ones obtained in Ref. [5] in the case of Nmax = 6 and Qq = 20 GeV^, taking into account 
a more approximate model for (n) and Fade approximations for CpJ (n) at n < 10. 

• At the scale Qq=20 GeV^ the obviously visible difference with the findings of Ref. [5] is related to the switch from 
the Fade approximant estimates of N'^LO contributions to coefficient functions to the expressions for Cp/'lint, 

/ON 

presented in Table 2, and obtained from the calculations of CpJ (n) at n = 1 [27] and n = 3, 5, 7, 9, 11, 13 [23]. 
The positive outcome of this change is the shift of from 340± 37 MeV to the the values, given in Table 6 ( 
which are less different from the results of the NNLO fits) and the minimization of the value of x^. 

• The LO and NLO fits seem to support the IRR model for the twist-4 terms by the foundation of negative values 

of A2 , which are different from zero within the presented error-bars and are in agreement with the results of the 
previous similar fits of Refs. [3, 5] and with the one, obtained in Ref. [38] using the NLO DGLAF analysis of the 
same set of CCFR'97 data. 

• At the NNLO the central value of A2 is also negative, but has large error bars. Moreover, the inclusion of 
the N^LO corrections clearly demonstrates the effective minimization of the free parameter of the IRR model, 
which becomes positive, but has statistical uncertainties twice as large as the central value. Thus, we may 
conclude, that at this level the interplay bewteen high-order perturbative corrections and the model for twist-4 
contributions, discussed from various points of view in Refs. [59]- [61], is manifesting itself. 



5.2 IRR approach and naive non-Abelianization 

Fcrturbativc expansion of the IRR model is usually understood within the framework of the large-/3o expansion, where 
/3o is the first coefficient of the QCD /3-function. The approximations for coefficient functions (but not anomalous 
dimensions) obtained within this model can be compared with the explicit expression for the quantities under con- 
sideration, calculated in the MS scheme. As a rule, the qualitative success of the estimating power of the large-/3o 
expansion is rather satisfactory (see e.g. [32, 53]). Let us study the application of this approach to the coefficient 
functions of odd Mellin moments of xFs used in our work. 

They can be presented in the following numerical form (see Refs. [27, 23]): 

c'-p^ = 1 -4^ -I- ^^(-73.333 -I- 5.333/) -h ^^(-2652.154 -I- 513.310/ - 11.358/^) (22) 

Cf^ = 1-M.667^-F A^(14.254- 6.742/) -I- ^^(-839.764- 45.099/ -M.748/^) 

Cf^ = 1 + 7.7A8A, + A^(173.001 - 19.398/) + A^(4341.081 - 961.276/ + 22.24:1 f) 

dp^ = 1 + 12.722^ + (345.991 - 30.523/) -h ^^(11119.001 - 1960.237/ + 43.104/^) 

&p^ = 1 + 16.915^ + ^2(520.006 - 40.355/) + (18771.996 - 2975.924/ + 63.171/^) 

C^^i) = 1 + 20.548^ + ^2(690.872 - 49.171/) + (26941.480 - 3984.412/ + 82.246/^) 

C^i^) = 1 + 23.762As -I- ^2(857.178 - 57.181/) + ^^(35426.829 - 4976.081/ + 100.351/^) 



Note that for the reasons discussed below, the order 0(A^)-corrections to Eq. (22) are presented without typical 
structures (wts), proportional to d"^'=d"'"^-terms. 

The procedure of large /Sg-expansion is formulated in our notations in the following way: one should extract from 
the explicit expressions for perturbative coefRcients, calculated in the MS-scheme, the leading terms in the number of 
flavours / and then make the substitution / — > — 6/3o, where /3o = /3o/4. This approximation is known in the literature 
as the "naive non-Abelianization" (NNA) procedure [62]. Note that it does not simulate (i"''^(i"''^-terms. Using this 
pattern we present below the coefficient functions of Eq. (22) in the NNA form: 





= 1 


-4^ + A2( 


-6 X 5.333/3o) + A^*(-36 x 11.358/3^) 




= 1 


+ 1.667^ + 


^2(6 X 6.742/3o) + 


^3(36 X 1.748/^2) 




= 1 


+ 7.748^ + 


Al{6 X 19.398/3o) H 


^Al{36 X 22.241/3^) 




= 1 


+ 12.722^ - 


h ^2(6 X 30.523/3o) 


+ Al{:m X 43.104/3^) 


^3 


= 1 


+ 16.915^ - 


1-^2(6 X 40.355/3o) 


+ A;^(36 * 63.171/3^) 




= 1 


+ 20.548^, - 


1- A^(6 X 49.171/3o) 


+ A^(36 X 82.246/3o) 


^3 


= 1 


+ 23.762^ - 


1-^2(6 X 57.181/3o) 


+ Al{36 X 100.351/3^) 



The obtained NNLO and N^LO corrections should be compared with the corresponding ones of Eqs. (24). The 

(2) (2) (2) (3) (3) (3) 

/-dependence of the ratios i?^?,^ ^ ya 

(n) = C';j{n)NNA/C';j{n) and R';^ 

NNAi^) — ^Fs {n)NNA/Cp^ (n)\u,is; which 
follow from the comparison of the related expressions of Eqs. (22,23), is presented in Table 7 and Table 8, where the 

/-dependence of the ratio Rp^{n)\]^Pade\ = C')?^ (n)|[i/i]/C)^^ (n)|u;ts is also given in round brackets. 



n 


/ = 3 


/-4 


/ = 5 


1 


1.25 


1.28 


1.31 


3 


-15.24 


-6.63 


-3.98 


5 


2.34 


2.61 


3.07 


7 


1.43 


1.42 


1.41 


9 


1.36 


1.41 


1.46 


11 


1.22 


1.24 


1.27 


13 


1.13 


1.14 


1.15 



Table 7. The /-dependence of the ratios Rpl nna{''^)- 



n 


/ = 3 


/ = 4 


/-5 


1 


1.61 (0.68) 


2.01 (0.87) 


2.99 (1.48) 


3 


-0.32 (-0.002) 


-0.26 (-0.1) 


-0.21 (-0.22) 


5 


2.46 (1.03) 


4.15 (1.38) 


41.32 (8.22) 


7 


1.4 (0.90) 


1.7 (0.99) 


2.4 (1.22) 


9 


1.1 (0.90) 


1.25 (0.96) 


1.56 (1.09) 


11 


0.95 (0.91) 


1.04 (0.96) 


1.20 (1.06) 


13 


0.85 (0.92) 


0.95 (0.97) 


1.02 (1.05) 



Table 8. The /-dependence of the ratios NNAi''^) Rp^ {n)\[Pade] (in round brackets). 

One can see that the NNA approach is working reasonably well in the case of the Gross-Llewellyn Smith sum rule 
(n = 1 moment). This fact was already observed in the review of Ref. [32]. It also gives satisfactory estimates both 
at the NNLO and N^LO in the case of odd values of n with n > 7, but does not work for n = 3 (where even the 
wrong sign is obtained) and for n = 5, where at the N^LO the subleading in / (and thus /3o) term is larger than the 
leading /3Q-contribution. Notice that in the case of even NS moments of F2 the situation was the same: the NNA 
approximation was predicting the correct sign starting from the n = 4 moment and was giving qualitatively good 
estimates in the cases of n = 6,8 moments [34]. Armed by the new information about explicit behaviour of the NS 
moments for F2 with n 10 [20] and n ~ 12, 14 [23], we extend the considerations of Ref. [34] to the case of higher 
moments, omitting / yZj^i e f-contribution to the order 0(>l^)-corrections of the NS moments of F2. Taking into 



account this approximation we present the explicit expressions for the coefScient functions of the NS moments of F2 
in the following numerical form: 

&plj^g = 1 + 0.444^^ + ^2(^7 694 - 5.333/) + ^^(442.741 - 165.197/ + 6.030/2) ^^A) 

&pIj^s = 1 + 6.607^ + ^2(141.344 _ 16.988/) + (4169.268 - 901.235/ + 23.355/^) 

dplj^g = 1 + 11.177^ + ^2(302.399 - 28.013/) + A^(10069.631 - 1816.323/ + 42.663/^) 

dpl^g = 1 + 15.530^ + ^2(470.807 - 37.925/) + A^(17162.372 - 2787.298/ + 61.9118/^) 

C^^°]^_5 = 1 + 19.301^ + ^2(639.211 - 46.861/) + ^3(24953.135 - 3770.102/ + 80.5201/^) 

dp^^j^s = 1 + 22.628As + ^2(804.585 - 54.994/) + ^^(33171.455 - 4746.441/ + 98.348/^) 

&p^^j^g = l + 25.611^s + ^2(9g5 _ 52.465/) + ^3(41657.116 -5708.216/ + 115.392/2) . 

The NNA versions of the expressions from Eq. (24) read: 

&pIj^s = l + 0.444^ + ^2(6 x 5.333/3o) + A3(36 x 6.030/3g) (25) 

&p1^s = 1 + 6.607^ + ^2(6 X 16.988/3o) + ^^(36 x 23.355/3^) 

&pIj^s = 1 + 11.177^ + ^2(6 X 28.013/3o) + A2(36 X 42.663/3^) 

Cflj^g = 1 + 15.530As + ^2(6 X 37.925/3o) + A^(36 X 61.9118/3^) 

Cj,"']vs = 1 + 19.301^ + ^2(6 X 46.861/3o) + A^(36 X 80. 5201/3^) 

^F^^NS = 1 + 22. 628As + yl2 (6 X 54.994/3o) + A;^ (36 X 98.348/3^) 

&p^^j^g = 1 + 25.611^ + ^2(6 X 62.465/3o) + ^1^(36 x 115.392/3^) 

The numerical values of the ratios of the coefficients of Eqs. (24,25), namely Rp^ nna{''^) — ^^f^ ns(''^)^^^/^f^ Nsi''^) 

and RF2,NNAi^) = CpJj^s{n)NNA/C^F2 Nsi'^)\wts are given below. In Tabic 10 Rp^i^^Ai^) is compared with 
Rp- {''^)\[Pade] = C")^ (n)|[]^/]^]/C)^, (r?,)|u,/,,. Getting support from the related results for n = 10. 12, 14 NS moments of 
F2 , we can make the conclusion that the findings of Ref . [34] and the new numbers for the moments of xFg, (see Tables 
7,8) demonstrate that at the NNLO and N^LO the NNA approximation is working in the NS channel for n = 1 and 
n > 6, which corresponds to the region of x, closer to the limit x = 1. In the cases of low NS moments the reason 
for failure of the NNA approximation remains unclear. A similar conclusion can be drawn from the comparison of the 
[1/1] Pade estimates for C^^-terms with the results of the interpolation procedure (see Table 2 and related discussions 

after it) and the /-dependence of the ratios Rp^ NNAi"^) Rps {''^)\[Pade] and Rp^ NNAi"^) Rf2 i'^)\[Pade] (see Tables 
8,10). 



n 


/ = 3 


/ = 4 


/ = 5 


2 


43.24 


-18.17 


-3.27 


4 


2.54 


2.89 


3.43 


6 


1.74 


1.83 


1.98 


8 


1.43 


1.48 


1.55 


10 


1.28 


1.29 


1.33 


12 


1.16 


1.18 


1.19 


14 


1.08 


1.09 


1.10 



Table 9. The /-dependence of the ratios Rp^ NNAi''^)- 
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-l.iO (-0.24j 
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2 54 fn 73"! 
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6 


1.55 (0.85) 


1.19 (0.93) 


2.74 (1.95) 


8 


1.18 (0.88) 


1.38 (0.94) 


1.71 (1.07) 


10 


1.00 (0.90) 


1.13 (0.95) 


1.31 (1.05) 


12 


0.90 (0.91) 


1.02 (0.96) 


1.09 (1.04) 


14 


0.82 (0.93) 


0.87 (0.97) 


0.95 (1.04) 



Table 10. The /-dependence of the ratios Rp^ nna{''^) -^^3 (?^)|[Pade]- 

This interesting similarity might be explained by the studies of the large- /3o hmit of the Pade approximant approach 
Ref. [63] and its relation to the BLM approach [64]. Note in these circumstances that within the large- /3o limit, the 
BLM approach was extended to all orders in perturbation theory in Ref. [65], while the possibility of incorporation 
of the subleading terms in number of flavours / and constructing the NNLO generalization of the BLM approach 
was first demonstrated in Ref. [66] (for further related analysis sec Ref. [67]). On the other hand the qualitative 
sucess of application of [1/1] Pade approximation for n > 6 might be considered as the additional argument in favour 
of possibility of its application for estimating N^LO contribution r(n) to the expanded anomalous dimension term 
Eq. (10). 



5.3 The determination of q.s{Mz) values and 
their scale-dependence uncertainties 



As is known from the work of Refs. [68, 35] it is rather instructive to consider the sensitivity of the results of the 

pcrturbativc QCD analysis to the variation of rcnormalization and factorization scales. We will study the question of 
factorization-renormalization scale dependence within the class of MS-like schemes. This means that we will change 
only the scales without varying the scheme-dependent coeflacients of anomalous dimensions and /3-function. 



The arbitrary factorization scale enters in the following equation: 



A(qVm|^) = a,(q7m^) 

where ji'^p is the factorization scale and 



1 + fci^(Q7M^) + fc2^j(Q7MF) + fc3A^(Q7MF) 



2 

MS ■■ 



Po 



5/3; 



kt + ^^kl 



/3o 



ki 



Let us choose the factorization scale as 
Then we have: 



MS 



ki = —Poln{kF) 



(26) 



(27) 



(28) 



In this case after application of the rcnormalization group equation and substitution of Eq. (26) into Eqs. (9,10) 
of Sec. 2 we get 



exp 



-dx 



= {A,{Q'/kF)r X AD{n,A,{Q'/kF)) 



where a = /2/3o and 



AD{n,A,{Q^/kF)) = 1 + 



p{n) + aki 



As{Q^/kF) 



(29) 



(30) 



+ 



q{n)+p{n)ki{a + 1) + ^kia + ^ ^ ' kl 
Po 2 



A'AQVkF) 



+ 



r{n)+q{n)ki{a + 2)+p{n) (§-A;i(a + l) + 

VPo 



(a+l)(a + 2) 2 



/3o 



kia+ ^klaQ + a) + 
Po ^ 



a{a+l){a + 2) 



kl 



Now let us consider the factorization and renormalization scale dependence, fixing ku = kp = k. In this case we 
should also modify the coefficient function in Eq.(6) as 



l + {n)A,{Q^/k) 



C]^>{n)-C'^;>{n)Poln{k) 



AliQ'/k) 



(31) 



+ 



4? W + 4? W iPoln^k) - (3iln{k)) - 2Cfjoln{k) 



AiiQ'lk) 



The commonly accepted practice is to vary k in the interval 1/4 < fc < 4 (see, e.g., Ref.[68]). We repeated our 

fits both without and with the IRR model of the twist-4 terms in the cases of A; = 1/4 and fc = 4. As in the fits 
described above, in order to achieve the minimum in we supplemented the interpolation procedure of the NNLO 
approximation for ^^p^ {A^) by fine-tuning of even terms 7j^^(6), (8) and 7^^^ (10) and got their values, comparable 
within small error-bars with the numbers given in Sec. 4. The same procedure was used in the process of the N'^LO fits. 
In fact they have more theoretical uncertainties than the NNLO ones. Indeed, in this case we applied the interpolation 
procedure to determine not only the NNLO coefficients of anomalous dimensions of even moments of xFt, SF, but the 
related N'^LO terms of the coefficient functions as well (see Table 2). The N'^LO corrections r{n) to AD-function in 
Eq. (30) were modelled using the [1/1] Fade approximant procedure. Note that for even values of n the numerical 
expresssions for q{n), which enter into the Fade approximants, are determined in part by the NNLO coefficients of 
anomalous dimensions of even moments of xF-^. In the process of N^LO fits in the case of fc = 4 the ambiguities 
of the applications of the [1/1] Fade approximation procedure reflect themselves in the necessity of supplementing 
the interpolation procedure by fine-tuning of the coefficient ^p' [12) in addition to the n = 6, 8, 10 NNLO anomalous 
dimension terms. Only after this additional step were we able to achieve a reasonable value of in this case also. 

The consequences of the study of factorization/renormalization scale dependence at the NLO, NNLO and N^LO 
in the case of the initial scale Qq = 20 GeV^ are presented in Table 11. 



Order 


N 

^ ^ max 


k 


All! 

MS 




A'2 {HT) 


/points 


NLO 


9 


1 


345±38 






85.1/86 




9 


1 


376±39 




- 0.121±0.052 


79.5/86 




9 


1/4 


482±57 


137 




90.0/86 




9 


1/4 


579±62 


203 


- 0.184±0.054 


78.8/86 




9 


4 


270±25 


-75 




84.7/86 




9 


4 


271±24 


-105 


- 0.032±0.051 


84.4/86 


NNLO 


9 


1 


331±37 






75.8/86 




9 


1 


332±36 




-0.015±0.051 


75.7/86 




9 


1/4 


379±45 


47 




78.7 /86 




9 


1/4 


399±46 


66 


- 0.084±0.052 


76.1/86 




9 


4 


297±27 


-35 




79.4/86 




9 


4 


318±30 


-15 


+ 0.117±0.052 


74.9/86 


N^LO 


9 


1 


327±34 






76.4/86 




9 


1 


329±34 




-f0.033±0.052 


76.0/86 




9 


1/4 


355±39 


28 




75.9 /86 




9 


1/4 


357±39 


28 


- 0.026±0.051 


75.9/86 




9 


4 


312±24 


-15 




74.8/86 




9 


4 


318±24 


-11 


+ 0.058±0.052 


84.5/86 



Table 11. The outcomes of NLO, NNLO and N^LO fits to CCFR'97 .TF3 data for > 5 GeV^ with different 
values of factorization/renormalization scales. The difference in the values of A^l^ is determined by lS.]^{MeV) = 

h_^{k) - h-^ik = 1). The value of the IRR model coefficient is given in GeV^. The initial scale is fixed at 
Ql = 20 GeV^. 



It should be mentioned that despite the approximate nature of the Pade resummation procedure used for the 
estimation of the N^LO contribution r(n) to the anomalous dimension function AD, in the case of application of the 
fine-tuning procedure at fc = 1 and k = 1/4 we get for jpj{ri) with n = 6,8, 10, 12 the numerical expressions which 
agree with the interpolated numbers of Tabic 1 in the 4th significant digit. For the fine-tuned fits with k = A and 

high-twist terms included the price for small values of is paid by more approximate determinations of 7^'' (10) 
(2) 

and jpj{12) which differ from the ones presented in Sec. 4 in the 3rd significant digit. For example, for n = 12 the 

l'2\ 

following number was obtained: 'ypj{12) w 1694.4907 ± 3.1194. It should be stressed, however, that this difference 

leads to negligibly small effects in the overall contributions to Pade estimated values of r{n). Moreover, it is rather 
pleasant that the approximate character of the fixation of this part of the theoretical input of our new N^LO fits does 
not drastically spoil reasonable (from our point of view) estimates for Jp^{n) at n = 6, 8, 10, 12. 

We make now several conclusions which come from the results presented in Table 11. 

• At NLO, NNLO and N^LO and fc = 1/4 the values of and thus are larger than in the case of A; = 1, 
while for /c = 4 smaller numbers are obtained. 

• The NLO and NNLO results of Table 11 are in satisfactory agreement with the similar ones from Table 6 of 
Ref. [5] , provided one takes into account the difference in the definitions of the parameter k (in Ref. [5] the case 
A; = 4 (fc = 1/4) corresponds to the choice A; = 1/4 (A; = 4) in Table 11). 

• To our point of view the results of the NNLO fits with A; = 1/4 (A; = 4) both without and with twist-4 terms 
simulate in part the results of the NLO (N^LO) fits with A; = 1. 

• The increase of order of perturbative theory approximations leads to minimization of the scale-dependence 
uncertainty which manifests itself through the decrease of the values of deviations. 

• In the case of NLO fits with HT terms the value of | | is larger than in the case of switching of power-supressed 
terms. However, this difference is minimized at the NNLO and the the N^LO especially. We think that this 
property is reflecting the correlation with the effective minimization of the fitted value of the HT parameter A2 , 
which becomes comparable with zero in the NNLO and N^LO fits with A; = 1. 

• We checked that for A; — 1 and A; = 1/4 the results are not sensitive to the changes of the initial scale from 

Ql = 20 GeV2 to Q2 = 5 GcV^. 

• However, when A; = 4, this pleasant feature is violated in the results of the NNLO and N^LO fits especially. 
Indeed, these fits are accompanied by the increase of up to over the 100/86 level. This fact can be related 
to pushing the value of Qq out of the considered kinematical region > 5 GeV^ (at Q§ = 5 GeV^ and A; = 4 
the region of 1.25 GeV^ < < 5 GeV^ should be also taken into account; however the NNLO and N^LO 
corrections to the coefficient functions are rather large in this region) . 

• It should be stressed that the NLO and NNLO results of Table 11 with A = 4 are closer to the ones obtained in 
the recent work of Ref. [37]. We will comment on the possible consequences of this observation in Sec. 6. 2 below. 

Let us now turn to determination of the values of as{Mz) from the results of Table 11. We transform into 
using the NLO, NNLO and N^LO variants of equation of Ref. [69]: 
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(35) 



where /?/ {p(^^) are the coefRcients of the /3-function with / (/ + 1) numbers of active flavours, Lh 
and Af/+i is the threshold of the production of a quark of the (/ + 1) flavour and C3 = -(80507/27648)C(3) 
(2/3)C(2)((l/3)Zn2 + 1) - 58933/124416+ (//9)[C(2) + 2479/3456]. 



These formulae contain the NNLO correction to the matching condition with coefficient C2 = —7/24, previously 
derived in Ref. [70] and correctly calculated in Ref. [71]. In our massless analysis we will take / = 4 and rrn, « 4.8 GeV 
and vary the threshold of the production of the fifth flavour from M| w m\ to M5 w {'orribf' in accordance with 
the proposal of Ref. [72]. The latter choice is based on the calculations of the LO and NLO massive corrections to 
the Gross-Llewellyn Smith sum rule. The final values of as{Mz) will be fixed at the middle of the interval, limited 
by the choices of threshold matching point at M5 w m\ and M| w (6mb)^. The appearing theoretical ambiguities 
reflect the uncertainties due to the manifestation of the massive-dependent contributions to the moments of xFy, in 
the massless fits. Another procedure of fixing the massive-dependent ambiguities in as{Mz), which result from the 
fits to xF^ data, was proposed in Ref. [73]. It is based on the application of the massive-dependence of the LO 
contribution to /3-fimction in the MOM scheme, previously studied in the works of Refs. [74, 75]. This procedure 
gives the estimates of the infiuence of mass-dependence on the value of as{Mz), extracted from CCFR'97 xFj, data, 
which are comparable to ours. To be more complete at this point, we also mention several other works, which are 
dealing with different prescriptions for estimating threshold uncertainties (see Refs. [76, 77] and the work of Ref. [78] 
especially, where massive-dependence of the MOM-scheme coupling constant was evaluated at the 2-loop level). It was 
shown in Ref. [72] that the application of the MS-scheme matching condition with the matching point M5 w (6mb)^ 
does not contradict the application of the massive dependent approach of Ref. [78]. Therefore, we can conclude that 
our estimates of massive-dependent uncertainties in as{Mz) can be substantiated by this comparisons of the results 
of Refs. [72, 78]. 

Taking into account the numbers given in Table 11, which were obtained with the twist-4 contribution fixed using 
the IRR model of Ref. [29], and the theoretical expressions of Eqs. (14)-(16) and Eqs. (32)-(35), supplemented with the 
estimates of the uncertainties due to different possibilities of the choice of matching point and experimental systematic 
errors, which come from separate consideration of this type of experimental uncertainties of CCFR'97 data, we arrive 
at the following values of as{Mz)-, extracted from the fits to CCFR'97 data for xF^ performed in this work: 



NLO HT of i?e/.[29] as{Mz) = 0.120 ± 0.0022(sfai) ± 0.005(sysi) 



(36) 



±Q.m2{thresh.)t°d^o^{scale) 



NNLO HT of Ref. [29] as{Mz) = 0.119 ± 0.002(stat) ± 0.005(syst) (37) 

±0m2{thresh.)toZ2i.scale) 

Minor differences with the similar results of Ref. [5] are explained by the incorporation of more significant digits in 

the process of calculations and by more careful study of the scale-dependence uncertainties. These values presented 
in Eqs. (36,37) should be compared with the one given by our new N^LO approximate fit: 

N^LO HT of i?e/.[29] as{Mz) = 0.119 ± 0.002(stai) ± 0.005(s2/si) (38) 

±Qm2{thresh.)tQ^m{scale) . 

Note again, that the experimental systematic uncertainties are extracted from the CCFR'97 data and were not taken 
into account in the process of our concrete studies for the reasons discussed above. As to the theoretical uncertainties 
of the as(Mx)-value, the incorporation of the high-order corrections to the fits leaves threshold ambiguities at the 



same level, but decreases scale-dependence uncertainties drastically ^. 

To study the influence of the twist-4 contributions of Eq. (21) to the values of as{Mz) we also extracted from 
Table 11 the corresponding results, obtained from the twist-4 independent fits to CCFR'97 data: 

NLOas{Mz) = 0.118 ±0.002(sfai)± 0.005(s2/s<) 
±0.002(i/ires/i.)t^:oo^(scak) 

NNLOas{Mz) = 0.119 ± 0.002(stat) ± 0.005(syst) (40) 
±0.002(t/ires/i.)t°;°|]^(.scale) . 

The approximate N'^LO twist-independent fits give us the following numbers: 

N^LOas{Mz) = 0.119 ± 0.002(.siaf) ±0.005(s2/si) (41) 
±Qm2{thresh.)to"om{scale) 

Let us make several conclusions, which follow from the comparison of the results of Eqs. (36)-(38) and Eqs. (39)-(40). 

• In the case when HT-corrections are included, the general tendency (as (M2))jvLO > {o-s{Mz))nnlo > {oisiMzYjmLO 
for the central values of the outcomes of the fits takes place. 

• The scale dependence of the NLO and NNLO results with HT-corrections included are larger than in the case 
of as(M2)-values, obtained without HT-terms. 

• The scale-dependence of the results of the approximate N^LO fits both without and with HT-corrections is 
almost the same. This feature is related to effective minimization of the contributions of HT-terms at the N^LO. 

• Starting from the NNLO the systematical experimental uncertainties dominate theoretical ambiguities as esti- 
mated by us. 

• The uncertainties of the matching conditions dominate scale-dependence ambiguities at the N^LO only. This 
might mean that the approximation of massless quarks works reasonably well in the analysis of CCFR'97 xF^ 
data up to NNLO. 

It is worth making several comments on the comparison of our results for as{Mz), extracted from the CCFR'97 
xF^ data, with those available in the literature. Within existing theoretical error bars, which reflect in part the 
special features of the procedures of taking into account threshold effects, our twist-independent NLO result is in 
agreement with the value as{Mz)NLO ~ 0.122 ± 0.004, given by the independent NLO fits to CCFR'97 xF^ data 
using the Jacobi polynomial method [73]. It is also in agreement with the one presented by the CCFR collaboration, 
namely as{Mz)NLO = 0.119 ± 0.002(exp) ± 0.004(t/ieorj/) [7], which was obtained with the help of the DGLAP 
approach and reproduced in Ref. [38] without treating carefully theoretical uncertainties. Within experimental and 
theoretical error bars our results also do not contradict the recent application of the Bernstein polynomial method 
[36] for the extraction of as{Mz) from CCFR'97 xF^, data at the NLO and NNLO [37]. Indeed, at the NLO it gives 
as{Mz)NLO = 0.116 ± 0.004, while at the NNLO they got as{Mz) = 0.1153 ± 0.0041(exp) ± 0.0061 (tfteor) [37]. It is 
worth mentioning here that despite the qualitative agreement with our results, the central values of as{Mz) obtained 
in Ref. [37] are lower than the central values of all existing NLO and NNLO determinations of as(Mz) from CCFR'97 
xF^ data. In Section 6 we will present more detailed comparison of the results of Ref. [37] with the ones obtained in 
our work and will propose a possible explanation of the origin of these deviations. 

Although we do not include in our fits the simultaneous analysis of the statistical and systematic experimental 
uncertainties, we think that our analysis has some advantages over other fits to CCFR'97 xF^ data. Indeed, the 
results of Refs. [73, 37] are free from considerations of the scale-dependence ambiguities studied in our work, while the 
DGLAP fits of the same data, performed in Refs. [7, 38, 39], do not take into consideration the contributions of the 
NNLO and N^LO perturbative QCD corrections, which we were able to treat in the way described above. 

To our point of view, another advantage of our work is that using more rigorous information than previously 
(see Ref. [5]) on the effects of the higher-order perturbative QCD contributions to the characteristics of xF^ SF, we 
continued the studies of the influence of twist-4 corrections to the extraction of as{Mz) from CCFR'97 xFs data. 

However, the IRR approach used by us is not the only way of modelling twist-4 effects. In the next Section we are 
considering the case when the twist-4 contribution will be approximated in a less model-dependent way. 

■^However, like in other studies of the CCFR'97 data of Refs. [3]- [5], [7], [37]- [39], [73], we arc neglecting theoretical uncertainties, which 
arise from the separation from the CCFR'97 data heavy nuclei corrections. Definite theoretical considerations of this problem [79, 80] were 
incorporated in the fits in Ref. [81] and indicate the decrease of as(Mz) at the NNLO level to the amount about 2 x 10~^. 



(39) 



5.4 Determination of the x-shape of the twist-4 corrections 



We now turn to a pure phenomenological extraction of the twist-4 contribution h{x) to Eq. (18). In order to get its 
x-dependence we will model h{x) by free parameters hi = h{xi), where Xi are the points in experimental data binning. 
The results are presented in Table 12. 





LO 


NLO 


NNLO 


N«LO 


X'^/nep 


69.1/86 


68.5/86 


65.6/86 


66.1/86 


A 


4.32 ± 1.34 


4.08 ± 1.13 


5.06 ± 0.46 


5.55 ± 1.27 


b 


0.629 ± 0.096 


0.616 ± 0.084 


0.682 ± 0.030 


0.711 ± 0.079 


c 


4.28 ± 0.14 


4.16 ± 0.15 


3.88 ± 0.20 


3.73 ± 0.34 


7 


1.91 ± 1.20 


1.87 ± 1.09 


0.73 ± 0.30 


0.25 ± 0.81 


Ag^ [MeV] 
aeiMz) 


327 ± 149 


395 ± 151 


391 ± 159 


370 ± 131 


U.i4U_o.oiO 


r. 1 90+0.007 
U.lZO_o 007 


u.iZ4_o 010 


r. -1 90+0.007 
U.lZO_o 008 




h{xi) 


GeV'^] 


0.0125 


0.016 ± 0.295 


0.056 ± 0.281 


0.054 ± 0.274 


0.072 ± 0.271 


0.025 


-0.055 ± 0.241 


0.038 ± 0.235 


0.239 ± 0.203 


0.284 ± 0.228 


0.050 


-0.079 ± 0.161 


0.113 ± 0.206 


0.425 ± 0.275 


0.472 ± 0.306 


0.090 


-0.231 ± 0.112 


0.050 ± 0.193 


0.146 ± 0.283 


0.158 ± 0.277 


0.140 


-0.369 ± 0.108 


-0.078 ± 0.132 


-0.012 ± 0.211 


0.013 ± 0.223 


0.225 


-0.492 ± 0.223 


-0.281 ± 0.160 


-0.038 ± 0.123 


0.062 ± 0.163 


0.350 


-0.344 ± 0.371 


-0.347 ± 0.367 


-0.207 ± 0.303 


-0.064 ± 0.227 


0.550 


0.129 ± 0.304 


-0.026 ± 0.335 


-0.172 ± 0.385 


-0.140 ± 0.325 


0.650 


0.398 ± 0.195 


0.275 ± 0.219 


0.144 ± 0.282 


0.131 ± 0.255 


n 








7F,^(n) 


6 






1247.9 ± 0.7 


1247.7 ± 0.6 


8 






1419.9 ± 0.8 


1419.9 ± 0.8 


10 






1561.9 ± 2.2 


1561.6 ± 2.2 


12 






1676.7 ± 5.3 


1676.8 ± 4.9 



Table 12. The values for of the the parameters h{xi), A, b, c, 7 and ALI with corresponding statistical errors. They 
are obtained from the fits with N^ax = 9 and Qq = 20 GeV . 

The x-shapes of h{x), obtained at LO, NLO, NNLO and approximate N^LO, are depicted at Fig.l, where we also 

illustrate the similar behaviour of h{x), obtained in the cases of N^ax = 6, 7, 8, which correspond to smaller number 
of Mellin moments, used in the perturbative part of the Jacobi polynomial reconstruction formula of Eq. (18). 

It should be noted that to minimize correlations between the values of hi, the parameters of the model for 
xF3{x,Ql) = A(Qg)x''('3o)(l - a;)=('3o)(l + -f{Ql)x) and the QCD scale A^, we choose 9 twist-4 parameters hi 
only, contrary to 16 ones considered in the process of our previous analysis of Refs. [3, 5]. The results of the fits are 
presented in Table 12. The approximate N^LO fits are based on the application of available N^LO corrections to the 
coefficient functions of odd Mellin moments for xF^ [23], supplemented with the smooth interpolation procedure (see 
Table 2) and using the [1/1] Fade model of the N^LO contributions r(n) to the expanded anomalous dimension term 
of Eq. (10) (see Table 3). In the process of NNLO and approximate N'^LO fits we faced a problem identical to the one 
revealed while fitting xFs CCFR'97 data without twist-4 terms (see Sec. 4) and with twist-4 contributions, modelled 
by means of the IRR approach (see Eqs. (20,21)). Indeed, to get the stable value of at Nmax = 9 it was necessary 

(n) 

to apply the fine-tuning procedure for the NNLO corrections to at n = 6,8, 10, 12. The obtained values of these 
parameters are presented in Table 12 also. Within error-bars they are in agreement with the numbers fixed by the 
smooth interpolation procedure (see Table 1). 

Several comments are now in order. 

1. The obtained values of x^, given in Table 12, are considerably smaller than the ones obtained in the process of 
the fits without twist-4 contributions (see Table 4) and with l/Q^-terms modelled through the IRR approach 
(see Table 6). This is the welcome feature of the analysis of DIS data, which is based on the model-independent 
parametrization of the twist-4 terms. 
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Figure 1: The x-shapcs of h{x) extracted from the fits to xF^ CCFR'97 data in different orders of perturbative theory. 
The initial scale is chosen at Qq = 20 GeV^. The cases of different Nmax are considered, a) Nmax=^'i b) A^maa:=7; c) 



2. The parameters A, b, c, 7 and given in Table 12 are in agreement with their values, obtained in Ref. [5] 
with the help of the fits, which were made in the case of N^ax = 6 and 16 HT parameters hi. 

3. Due to the effect of correlations of hi and the values of the QCD scale parameter have rather large statistical 
error-bars. 

4. At the LO and NLO the a;-shape for h{x) is rather stable to the increase of N^ax from N^ax = 6 to N^ax = 9 
and therefore, to the incorporation of the additional Mellin moments in the procedure of reconstruction of xF^ 

via the Jacobi polynomial technique. 

5. The x-shape of h{x), obtained at the LO and NLO, is in agreement with the prediction of the IRR model of 
Ref. [29]. 

6. At the NNLO we observe the sinusoidal-type oscillations of h(x), which are becomimg more vivid in the cases 

of Nmax = 8 and N^ax = 9. 

7. This feature seems to be in agreement with the qualitative expectations which result from an educated guess 
about the possible modifications of the prediction of the IRR model at the NNLO [82]. 

8. At the N^LO, for all considered Nmax, the shapes of h{x) are similar to the ones obtained in the process of new 
NNLO fits. 

9. It is worth mentioning that the positive bumps in the x < 0.1-region of the NNLO plots of Fig.l appear after 

(n) 

applying the fine-tuning procedure to NNLO coefficients of ^ for n = 6, 8, 10, 12, which gives us the possibility 



to get reasonable values for both and ALl at NNLO and beyond. In otherer words, the x-profile of the twist-4 



contribution is related to the values of the NNLO corrections to 7^^'' for even n. We consider this observation 
as an additional argument in favour of getting explicit results for these terms. 



10. In general, taking into account systematic experimental uncertainties of the CCFR'97 data for xF^ might make 
the sinusoidal- type oscillations of h{x), which demonstrate themselves in the NNLO and approximate N^LO 
fits with smaller number of free parameters hi, less vivid and more comparable with zero. Indeed, the effective 
minimization observed previously in Refs. [3, 5] of the contribution of twist-4 terms during less definite than 
present NNLO fits to xF^ CCFR'97 data were confirmed in the process of the NNLO DGLAP fits to F2 data 
for charged leptons DIS [83, 84], which were based on the application of the approximate NNLO model for the 
DGLAP kernel from Ref. [68]. The observed changes in the a;-shapc of twist-4 contributions h{x) to xF^ and F2 
SFs serve as additional arguments in favour of high- twist duality, which demonstrates itself through the interplay 
between NNLO perturbativc QCD corrections and 1/Q^ terms. 

11. It is interesting to note, that large assosiated with considerable twist-4 contributions was revealed some time 
ago in the process of NLO DGLAP fits of other vN DIS data [85]. Note, however, that in this work twist-4 
contributions were not free, as in our analysis, but simulated using special model. 



6 Comparison with the results of other NNLO analyses 

In this Section we compare the results of our studies with the outcomes of other analyses of arFs at the NLO and 
beyond, pcformcd independently in the works of Ref. [35] and Ref. [37]. 

In Ref. [35], to study the evolution of the NS contributions to F2 and xF^ up to the approximate N"^LO level 
of massless QCD, the NNLO corrections to the DGLAP equation coefficient functions for xF^, [13] and the N^LO 
corrections to definite Mellin moments of NS SFs, obtained in Refs. [20, 23], were combined with the NNLO model 
for the NS kernel, previously obtained in Ref. [68]. 



6.1 Comments on estimates of scale-scheme dependence 
uncertainties. 



Using the input 



F2,Ns{x,Ql) = xFsix, Ql) = x°-\l - xf 



(42) 



specifying the reference scale Qq through the normahzation condition as{ii%^ = Qq = 30 GeV^) = 0.2 , irrespective of 
the order of the expansion and varying the renormaUzation scale in the conventional interval jQ^ < n\ < 4(5^, the 
authors of Ref. [35] studied the effects of the scaling violation using a definite model for the xF^^ data. The values 
of as (30 GeV^) for / = 4 are given in the second column of Table 13. However, while analyzing real CCFR'97 xF^ 
data, we found that scale-dependent uncertainties can be larger. Indeed, using those numbers from Table 11, which 
are related to twist-independent fits to the CCFR'97 data, we obtain the inputs in the third column of Table 13: 



Order 


Values from Ref. [35] 


Our values 


NLO 
NNLO 
N^LO 


0.2035tl!;lil^ 
0.1995t°;°°?^ 
0.2000t°;°°g^ 


0.2104+--5^ 
0.2148t°:rr^ 
0.2144troS^^ 



Table 13. The comparison of the scale-uncertainties of as(30 GeV^), obtained in Ref. [35] and in the process of our 
studies. 

At the qualitative level both sets of numbers are in agreement with each other. Moreover, the scale-dependence 

of both sets of numbers has the tendency to decrease from NLO up to N^LO. Definite differences between the central 
values of as may be traced to the fact that our results for as(30 GeV^) correspond to as{Mz) ~ 0.118 irrespective 
of the order of the expansion (see Eqs. (39)-(41)), while the choice as(30 GeV^) = 0.2 in Ref. [35] corresponds to a 
lower value as{Mz) » 0.116. 

On the other hand the scale-dependence of our twist-independent NLO and NNLO results for as{Mz) (see 
Eqs. (39), (40)) is in agreement with the previous estimates of this kind of theoretical uncertainties, namely 

^as{Mz)NLO=tl-Zl > AQ«(Mz)iviVLO =iSi (43) 
obtained in Ref. [68] using the model constructed in this work for the NNLO NS DGLAP kernel. 

As to the application in Ref. [68] of the renormalization-scheme optimization methods of Refs. [47, 86] for estimating 
higher-order corrections (up to N^LO) to the factorization-scheme independent quantity , defined as 

dlnM^3 aC^"V^)/aA, 

we think that it might give larger theoretical uncertainties than those presented in Ref. [68]. Indeed, we previously 
used this ratio in Ref. [5] in the process of the attempt to perform the massless NNLO fits to CCFR'97 xF^ data using 
the effective-charges (ECH) approach of Ref. [87] (for the related methods see Refs. [88]- [91] and the independent 
unpublished proposal of Ref. [92] ; for the related phenomenological applications in the NLO fits to the charged leptons 
DIS SFs data see Ref. [93]). As was found in Ref. [5], the NNLO ECH fits to CCFR'97 xF^ data face the problem 
of the drastical increase of up to the level of ~ 111/86. This effect was explained by the appearance of large 
and negative values of the NNLO coefficients jS^^'' of the ECH /3-functions, related to the scheme-invariant quantities 
introduced within the context of the principle of minimal sensitivity (PMS) Ref. [94] . Here we are going to demonstrate 
explicitly how this effect appears, using the results presented in Sec. 3. 

In the case of n > 2 the NNLO approximation of the factorization-scheme independent kernel Kn can be rewritten 
in the following form: 

Rn = =As + d, {n)Al + d2 {n)Al (45) 

Putting di (n) = d2 (n) = we arrive at the following renormalization-group equation 



l3i}\{Rn) = -2(/3oRl + /3,Rl + pi''^Rt) , (47) 



with the effective /J-function defined as 

dRn ^ 
^ dfj, 

where the NNLO coefficient of Eq. (47) is related to the NNLO coefficients /32 of the MS-scheme /^-function as 

/32+A(n) (48) 



with 

A{n) = Po(^d2{n) - n^in)^ (49) 

and 

n2{n)=di{n)(^^+di{n)j . (50) 

The similar equations can be derived at the N^LO and beyond (see Refs. [47, 86]). 

The ECH-inspired estimates proposed in Refs. [47, 86] work only in the case when the differences /S^"'' — Pk with 
k > 2 are small. These conditions turned out to be valid for the e"'"e~-annihilation Adler D-function, DIS sum rules 
in QCD [47] and {g — 2)^ in QED [86] also. But unfortunately, they are not working in the case of the quantity 
defined by Eq. (45). Indeed, using the results from Tables 1,2 we obtain the numerical expressions for A(n), which 
are presented in Table 14: 



n 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


A(n) 


-1976 


-1288 


-1066 


-937 


-851 


-783 


-730 


-684 


-644 


-606 


-576 


-547 



Table 14. The n-dependence of A{n) = /j^"^ - P2 for / = 4. 

Comparing now these numbers with the numerical value for (^2 for / = 4. (/?2 = 406.35), wo arrive at the conclusion 
that the basic assumption of the ECH-inspired estimates of Refs. [47, 86], namely P2^^ « /32, does not work at NNLO, 
not only in the case considered previously of the correlator of scalar quark currents [53], but for the factorization-scheme 
independent quantity of Eq. (46) as well. 

Despite the fact that we did not consider in detail the case of F2^ns' SF, we think that the estimates presented in 
Ref. [35] of perturbative theoretical uncertainties for the factorization-scheme independent kernel Kn for F2^ns made 
with the help of the ECH and PMS approaches at the NNLO and beyond might be underestimated (at least in the 
case of n < 13). However, the general tendency of the absolute value of A(n) to decrease with increasing n might lead 
to the improvement of the situation for a larger number of moments, which are limited in Ref. [35] by n = 30. This 
guess can be substantiated only after completing the explicit analytical calculations of NNLO corrections to DGLAP 
kernels. 

Another interesting subject, related to negative values of ^2"^ corresponds to the appearance of perturbative IRR 
zeros of the ECH /3-function at the NNLO, considered previously as spurious ones in a number of works (see Refs. [95]- 
[98], [53]). It should be stressed, that for the quantity Kn these perturbative IRR zeros are manifesting themselves 
obviously in the case of n = 2 and less obviously for n = 3. In the first case the critical value of the corresponding 
effective charge {as)eff is small, namely 0.4, while for n = 3 it is over 0.7, which is rather close to the non-perturbative 
region. For n> A these zeros lie in the typical nonperturbative sector where {as)eff > 1. In view of this it might still 
be possible to apply the ECH or PMS approaches for the analysis of the scheme-dependence of the NNLO perturbative 
QCD predictions for xF^ moments with n > 3. It could be of interest to study this problem in the future. 



6.2 Comments on outputs of the NNLO Bernstein polynomial 
analyses 

Let us now comment on the comparison of our NLO and NNLO results for A^ and as{Mz) with the ones obtained 
in another interesting work of Ref. [37]. The authors of this work used the theoretical input identical to ours, namely 
the results of the NNLO perturbative QCD calculations for the anomalous dimensions and coefficient functions of 
odd Melhn moments of xF^, and performed the NLO and NNLO fits to CCFR'97 data for xF^ with the help of 
the Bernstein polynomial technique, proposed in Ref. [36]. In the process of these fits, the initial parametrization 
xF^{Q'q) = Ax^{l — xY, which is similar to Eq. (6), was considered. The initial scales Qq = 8.75 GeV^ and Qq = 
12 GeV were chosen inside the kinematical region of CCFR'97 data 7.9 GeV^ < < 125.9 GeV^ which is only part 
of the region used in our work. 

The results of twist-independent fits at the scale Qq = 8.75 GeV^ obtained in Ref. [37] are summarized in the last 
two columns of Table 15, where the error bars include statistcial and systematic experimental uncertainties. 

These results can be compared with our outputs in Table 4, which result from our Jacobi polynomial twist-4 
independent fits of the CCFR'97 data for xF^ in the kinematical region Q2 > 5 GeV^ with the cuts W > 10 GeV^ 



X < 0.7. It should be stressed that this region is identical to the one studied in the original work of the CCFR 
collaboration [7] . In particular, we are interested in the following values of 

= 265 ± 36 MeV (51) 
= 347 ± 37 MeV 
= 332 ± 35 MeV 

which correspond to the choice of the initial scale Qq = 8 GeV^ (note that the results of Table 4 demonstrate that 
these values are almost independent of the choice of Qo)- 

It is quite understandable why our results have smaller error-bars: contrary to the results of Rcf. [37] they are 
defined by the statistical experimental errors of CCFR'97 data alone. However, the explanation of the discrepancies 
in the central values of A^—L- is not so obvious. 

MS 

In order to clarify the situation we performed the Jacobi polynomial fits for two sets of experimental data from 
the CCFR'97 collaboration, choosing the same initial scale Qq = 8.75 GeV^ as in Ref. [37]. 

1. First, we considered the same data set as in Ref. [37], i.e. the kinematical region 7.9 GeV^ < Q2 < 125.9 GeV^. 

• The comparison of the results of our twist-4 independent fits with the ones given in Table 3 of Ref. [37] are 
presented in Table 15 below. 



Order 


MS 


as{Mz) 


X^/nep 


Al^ [37] 

MS L J 


a^iMz) [37] 


LO 
NLO 
NNLO 


227±37 
298±38 
303±38 


0.1309±0.0037 
0.1169±0.0025 
0.1187±0.0026 


92/74 
76/74 
65/74 


217±78 
281±57 
255±55 


0.130±0.006 
0.116±0.004 
0.1153±0.0041 



Table 15. The values of A^ and as{Mz), obtained with the help of Jacobi and Bernstein polynomial 
techniques Ql = 8.75 GeV^ . 

Notice that the difference between the values of A^ obtained by the Jacobi and Bernstein polynomial 
techniques are minimized in this case. The LO and NLO results for as{Mz), as obtained by us, almost 
coincide with those taken from Ref. [37]. However, at the NNLO our value of as{Mz), which is in agreement 
with the result of Eq. (40) , is comparable with the similar one given in Ref. [37] only within the presented 

experimental error-bars, which in the latter case also include systematical experimental uncertainties. 

/ 1 \ 

• To perform a more detailed comparison we also estimated the uncertainties in the extraction of A^^^ at the 
NNLO, as those considered in Table 4 of Ref. [37]. The results are given in Table 16. 



Sourse of errors 


AilL 

MS 


AAlfL 

MS 


A^ [37] 

MS L J 


AAi^ [37] 

MS L J 


No TMC 


326 


23 


298 


43 


HT 


316 


13 


270 


15 


Ql to 12 GeV^ 


298 


-5 


263 


-8 


NNLO* 


294 


-9 


209 


-46 



Table 16. Theoretical uncertainties of A^, obtained with the help of Jacobi and Bernstein polynomial 
techniques Ql = 8.75 GeV^. 

In Table 16 we mark by the symbol NNLO* the uncertainties of A^IL obtained from the NNLO fits with 

as defined through its N^LO expression (see Eqs. (14)-(16)). One can see that we obtained twice as small 
uncertainties while neglecting TMC and five times smaller effects while introducing the N^LO expression 
for as in the NNLO fits. It should be stressed that the similar small difference between the the values of 
A^ of the NNLO and NNLO* fits to the CCFR'97 xF^ data was aheady observed in Ref. [5] for a larger 
kinematical region and for different values of Qq. 

• At present we are unable to explain the most significant differences with the NNLO results of Ref. [37] . We 
think that more detailed comparison of the Jacobi and Bernstein polynomials approaches at the NNLO is 
really on the agenda. 



NLO Ag) 
NNLO aW 



2. If, following the CCFR collaboration, we exclude one data point with < 10 GeV^ (namely the point with 
X = 0.65 and = 12.6 GeV^), which has a large systematical error, and if we take into account the complete 
data set with > 5 GeV^ and the cut W'^ > 10 GeV^, we reproduce the results of Eq. (51) with rather small 
theoretical uncertainties: 



Order 


AlfL 

MS 


as{Mz) 


X^/nep 


LO 

NLO 

NNLO 

NNLO no TMC 
NNLO* 


265±36 
340±37 
333±36 
360±32 
322±35 


U.iO*0_Q QQ33 
U.liaOQ QQ23 

1 206+0 008 

U.lZUO_g QQ20 

n 1 1 9Q+0.0018 

U.iiZO_g QQ20 

llQQ+0 0021 
u.±±aa_Q QQ22 


113/86 
87/86 
74/86 
77/86 
74/86 



Table 17. The results of the fits to the CCFR'97 data within the kinematical conditions, used in our work. 
Q§ = 8.75 GeV^ 

It is worth noting that theoretical uncertainties due to the omission of TMC and due to the consideration of the 
N^LO expression for as, Eqs. (14)-(16), in the NNLO fits remain the same, as in the case of Table 16. 

Several additional comments are now in order 

• In Table 15 and Table 17 the values of as{Mz) were obtained from Eqs. (32)-(34) using the matching point 
M5 = mb « 4.8 GeV. Thus we neglected the uncertainties due to fixation of 6-quark threshold ambiguities, 
described in Sec. 5.3. 

• Note, however, that our estimate {Aas{Mz))thresh ~ ±0.0017 is only slightly larger than the estimate 
{Aas{Mz))thresh ~ 0.0010, given in Ref. [37] after application of a different method for estimating 6-quark 
threshold uncertainties. 

• It should be noted that contrary to the analysis in Ref. [37] we were able to study the scale dependence 
uncertainties of our results in the case when the renormalization scale was taken equal to the factorization 
scale. The results of our studies, presented in Sec. 5.3 and Table 11 demonstrate the following interesting 
feature: in the case of A; = 4 and HT neglected, both NLO and NNLO results for are almost identical 
to the ones obtained in Ref. [37] in a narrower kinematical region of CCFR'97 xF^ data. Therefore, a possible 
explanation of the deviations of our results from the ones in Ref. [37] might be related to the fact that scale- 
dependence ambiguities of the latter were not studied and might increase the theoretical uncertainty for 
as{Mz). 

6.3 The model independence of high- twist duaUty effect 

Despite the fact that in Ref. [37] the subject related to the inclusion of twist-4 terms was briefly considered at the 
NNLO only, it is rather instructive to perform a similar analysis at the LO, NLO and repeat the NNLO studies using 
the .Jacobi polynomial approach. It should be mentioned that in Ref. [37] the more simple than IRR-model form of 
the twist-4 corrections was used, namely 

<J^3(Q') = n^^M^HQ') , with S; = a{A^f (52) 

Its coefficient function differs from C{n) of Ref. [29] (see Eq. (21)), which for the moments under consideration has the 
following numerical values: (7(2) = 1.6667, (7(3) = 1.6333, (7(4) = 1.4, (7(5) = 1.0381, (7(6) = 0.5857, (7(7) = 0.0659, 
(7(8) = -0.5063, (7(9) = -1.1205, C'(IO) = -1.7689 and (7(11) = -2.4461, etc. Notice that starting from n=8, C{n) 
changes the sign and in the asymptotic regime tends to — n. Therefore, it is of interest to investigate the model- 
dependence of the effect observed in Sec. 5.1 of high twist duality using the HT model of Eq. (53), which is different 
from the IRR model of Ref. [29]. 

This question was studied by us using the set of CCFR'97 xFs data considered in item (1) of Sec. 6. 2, The results, 
obtained with the help of Jacobi polynomial fits, are presented in Table 16, where 6-quark threshold uncertainties were 
not taken into account. Note that in our fits we considered i?2 ^■s the free parameter, and then determined the value 
of the parameter a, considered to be free in the fits of Ref. [37] . 



Order 


MS 


X^/ncp 


B'2 


a 


asiMz) 


LO 

NLO 

NNLO 


433±89 
371±72 
316±51 


88/74 
75/74 
64/74 


-0.330±0.126 
-0.135±0.113 
-0.031±0.088 


-L76±0.37 
-0.98±0.57 
-0.31±0.80 


+U.Uub6 
^0.0063 
n 1010+0. 0039 

U.1Z10_Q gQ44 

U.11»Oq QQ44 



Table 18. The results of the fits to the subset of CCCFR'97 data with HT contribution, considered in Ref. [37]. The 
initial scale is chosen as Qq = 8.75 GeV^. 

Comparing now Table 18 with Table 6 of Sec. 4 we arrive at the following observations: 

• In both cases the values of are almost the same and thus do not depend on the typical structure of the 
twist-4 model. 

• At LO and NLO the value of the parameter B2 is in agreement with the value of the parameter A2 of the IRR 
model and decreases after NLO effects are taken into account. 

• At NNLO the value of the parameter is a bit larger than the similar value of A2, but within error bars both 
are compatible with zero. 

• At NNLO the central value for the parameter a is a bit larger than its value obtained in Ref. [37], but within 

existing uncertainties they arc compatible. 

• The effect observed in Sec. 5.1 of interplay between perturbative QCD and twist-4 corrections remains also valid 
in the case of the choice of twist-4 model given by Eq. (52). 



7 Conclusions 



In this work we used the new perturbative QCD input in the form of NNLO corrections to the anomalous dimensions 
and N^LO corrections to the coefficient functions for odd xF^ moments [23] to improve our previous fits to xF^ 
CCFR'97 data, performed in the works of Refs. [3, 5]. We demonstrated that the application of the smooth interpolation 
procedure, supplemented with the fine-tuning of NNLO corrections to 7^^^^^ for several even n gives us the chance to 
include in the NNLO and approximate N^LO fits a greater number of Mellin moments up to n < 13. The basic feature 
which we revealed in the process of our fits is the drastical reduction of the scale-dependence uncertainties for as{Mz) 
at the NNLO and beyond. The obtained values of as{Mz) turned out to be in in agreement with the world average 
value of this parameter, namely as{Mz) ~ 0.118 [99, 100]. The previously discovered property of interplay between 
NNLO perturbative QCD corrections and twist-4 terms is confirmed using different models for the l/Q^ corrections. 
The feature observed in the process of the NNLO fits with model-independent parametrization of the twist-4 terms 
is the sinusoidal oscillation of its .x-shapc around zero with definite positive bumps in the low x region. At present 
we do not know whether this typical behaviour can be described by the NNLO generalization of the IRR model of 
Ref. [29], or whether it will disappear after taking into account systematic experimental uncertainties. This problem 
can be studied using the machinery of the work of Ref. [40], which can allow us to fix experimental uncertainties of 
as{Mz) and twist-4 terms on more solid ground. 

As to the effective decrease of the twist-4 contributions at the NNLO and beyond, we are unable to disfavour the 
possibility that it occurs because of the use of the CCFR'97 data for xF^, which are still not precise enough. Possible 
future more detailed DIS uN data, which are expected to be obtained at the Neutrino factory [101], might be useful 
for clarification whether it is reliable to detect more clear signals from twist-4 contributions at the NNLO level. 

On the other hand our NNLO Jacobi polynomial fits revealed the necessity of getting more precise (namely exact) 
values for the NNLO corrections to the anomalous dimensions for even n, which are related to still explicitly 
uncalculable NNLO corrections to the kernel of the DGLAP equation for xF^. Having this information at hand, one 
might be able to perform NNLO Jacobi polynomial fits avoiding interpolation and fine-tuning procedures, which were 

(n) 

used by us to fix NNLO corrections to 7]^^ for even n, and fix still remaining theoretical uncertainties in the x-shape 
of HT-contribution to xF^, extracted at the NNLO. More detailed understanding of other physical effects when using 
explicit NNLO corrections to the kernels of the DGLAP equations might be revealed in the process of more detailed 
comparisons between different methods, which implement the classical DGLAP solution of differential equations in 
the X space and the Jacobi and Bernstein polynomial techniques. 
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Erratum to the revised version of hep-ph/0106221 



Improved fits to the xF^ CCFR data at the next-to-next-to-leading order and beyond. 
By A. L. Kataev, G. Parente and A. V. Sidorov, 
pubhshed in 

Physics of Elementary Particles and Atomic Nuclei, vol. 34 (2003) pp. 43-87 
[Physics of Particles and Nuclei vol. 34 (2003) pp. 20-46] 

The bug crept into the calculations of the numerical values of the Ag-coefHcients Cp^ (n) in the QCD expression for 

the coefficient function C|,"^(A,) (the definit ion see in Eq.(8)) of the odd Mellin moments of Eq.(l) for xF^ structure 
function of deep-inelastic neutrino-nueleon scattering [Ag — ag/ {A/t:)) with n = 1, 3, 5, 7, 9, 11, 13 and f — A numbers 
of flavours. This bug resulted from using in the computer subroutine, which calculated the values for CpJ (n) from 
the given in Eq. (22) order 0{Al) approximations for Cp^\ where instead of in the last terms / was typed. These 

(3) 

errors affected also the values of the even CpJ (n) coefficients, obtained from the explicitly calculated ones using the 
smooth interpolation procedure. The corrected results are given below in the 5-th corrected column of Table 2 of the 
paper. 



n 












CK^WI[0/21 


1 


-4 


-52 


-52 


-780.6427 


-676 


480 


2 


-1.778 


-47.472 


(-46.4295) 


(-1206.83008) 


-1267.643 


174.4079 


3 


1.667 


-12.715 


-12.715 


-992.198975 


97.00418 


-47.01328 


4 


4.867 


37.117 


(37.0076) 


(-269.865143) 


283.0851 


246.0090 


5 


7.748 


95.4086 


95.4086 


851.838501 


1174.835 


1013.328 


6 


10.351 


158.2912 


(158.4032) 


(2286.68115) 


2420.569 


2167.903 


7 


12.722 


223.8978 


223.8978 


3967.71313 


3940.284 


3637.790 


8 


14.900 


290.8840 


(290.8421) 


(5844.3042) 


5678.657 


5360.371 


9 


16.915 


358.5874 


358.5874 


7879.04004 


7601.721 


7291.305 


10 


18.791 


426.4422 


(426.5512) 


(10044.4785) 


9677.391 


9391.308 


11 


20.544 


494.1881 


494.1881 


12319.7676 


11885.25 


11633.28 


12 


22.201 


561.5591 


(561.2668) 


(14687.1133) 


14204.22 


13991.80 


13 


23.762 


628.4539 


628.4539 


171728.1191 


16620.99 


16449.68 



Table 2. The values for NLO, NNLO, N'^LO QCD contributions to the coefficient functions, used in our fits, and 
the results of N^LO Pade estimates. 

The application of the corrected numbers in the next-to-next-to-next-to-leading order Jacobi polynomial fits of the 
experimental data of the CCFR-collaboration resulted on slight decrease of of N^LO values for as presented in 
Tables 6, 11 and 12 by 3 MeV only and does not affect any conclusions of the paper. 



